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Abstract 

We describe briefly new program for SUSY masses calculations. The main dis- 
tinction of our program is that we start to solve the renormalization group equations 
for soft SUSY breaking parameters for SU(5) SUSY GUT model from Planck scale 
Mpl = 2.4 • 10 18 Gev. Our program works also for large tan(/3). We find that for 
mo < 0.5 • mi the effects of the evolution from Planck scale to GUT scale are very 
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essential. In particular, we find that neutralino even for small rriQ is always LSP. We 
also introduced in our program some parameter of non-universality of the gaugino 
masses at GUT scale. Playing with the non-universality of the gaugino masses at 
GUT scale it is possible to have the situation when leptonic modes are suppressed 
and the single SUSY signature is hadronic jets with missing energy. 
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Supersymmetric electroweak models offer the simplest solution of the gauge hierarchy 
problem 0-0]. In real life supersymmetry has to be broken and the masses of superpar- 
ticles have to be lighter than 0(1) Tev [ffl. Supergravity gives natural explanation of the 
supersymmetry breaking, namely, an account of the supergravity breaking in hidden sec- 
tor leads to soft supersymmetry breaking in observable sector |3]]. For the supersymmetric 
extension of the Weinberg-Salam model soft supersymmetry breaking terms usually con- 
sist of the gaugino mass terms, squark and slepton mass terms with the same mass at 
GUT scale and trilinear soft scalar terms proportional to the superpotential M. 

In this paper we briefly describe our program "PLANCKSUSY" for SUSY masses 
calculations ( the program can be received by e mail: Vsevolod.Popov@cern.ch). As 
it has been mentioned in the abstract the main peculiarity of the program is that we 
can start not from GUT scale Mqut ~ 2 • 10 16 Gev but also from Planck scale Mpl = 
2.4 ■ 10 18 Gev. We assume that physics between Planck scale and GUT scale is described 
by standard supersymmetric SU(5) GUT model M. Of course, we don't know for sure 
the physics between Planck and GUT scales so we consider SU(5) SUSY model as the 
simplest possibility to use it as a " quasi- realistic" model to take into account the effects of 
the evolution between Planck and GUT scales. The standard assumption |3], H| is to start 
from GUT scale and to impose universal mass relations for squark, slepton and Higgs soft 
masses and for gauginos at GUT scale 

m sq = m s i = m H = m Q , (1) 

7713 = m,2 — mi — mi (2) 

2 

However in standard approach soft SUSY breaking terms result of the supergrav- 

ity breaking in hidden sector so it is more valid to impose universal boundary conditions 
not at GUT scale but at Planck scale. It should be noted that in superstring inspired 
models in general we expect non-universal boundary conditions for SUSY soft breaking 
terms To take into account possible effects from non-universal boundary conditions 
we introduced in our program some parameter of non-universality for gaugino masses at 
GUT scale. It should be noted that our program takes nonzero Yukawa couplings for 
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t-quark, b-quark and tau-lepton so we can use it for large values of tan(/3). Besides we 
can start both from Planck and GUT scales to compare the results of the calculations. 

Our program can be divided roughly speaking into three parts. At first stage we 
solve numerically renormalization group equations for SUSY soft breaking parameters 
and coupling constants for SU(5) SUSY model between Planck and GUT scales. Here we 
use one loop RG equations of ref . @ . At second stage we solve numerically RG equations 
for SU{3) ® 577(2) <g> U(l) MSSM between M GUT and M SUSY . We use two loop RG 
equations for couplings and gaugino masses and one loop RG equations for other soft 
SUSY breaking parameters 0. From Msusy to m t we use RG equations for standard 
WS model. At third stage we solve the equations for the determination of the electroweak 
minimum of the one loop effective potential to determine the absolute value of the (jl 
parameter and the Higgs boson masses. We also calculate mixing effects for neutralino, 
chargino and stops. We calculate the pole gluino mass and other running sparticle masses 
m(m). The difference between running and pole masses for squarks and other sparticles 
is rather small. Here we use the formulae of refs.|5|1-[Pl. The input parameters of our 
program are the standard ones (mo, mi, tan(/5). sign(fj,), A, mt). Our program can start 
both from Planck scale and GUT scale. We also introduced explicitly some parameter of 
the non-universality for gaugino masses at GUT scale. We use the following formula for 
gaugino mass matrix at GUT scale: 



where I is the unit matrix and $ = $oDiag(2, 2, 2, —3, —3) Parameter k determines the 
non- universality of gaugino masses at GUT scale. Our formulae for SU(3), SU(2) and 
U(l) gaugino masses can be written in the form 



m = mil + fc$ 



(3) 



777,3 = mi (1 + 28) 



(4) 



m 2 = mi (1 — 35), 
mi = mi (1 — 5) 



(5) 



(6) 
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Here 5 is some parameter of the gaugino mass non-universality. 
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Let us briefly consider first nontrivial physical consequences of our program. We have 
found that for m > 0.5 -mi both "Planck" and "GUT" spectra coincide up to 20 percent. 
For tjiq < 0.5 ■ mi we find essential difference in slepton spectrum. The most interesting 
difference is that for small m and mi > 0(300) Gev the lightest superparticle in "GUT" 
model is charged right-handed tau slepton that contradicts to the experimental data 
on abundances of anomalous super- heavy isotopes and usually this region in (m ,mi 
plane is considered to be "theoretically excluded". When we start from Planck scale 
the lightest superparticle in this region of (m , mi) parameters is always neutralino so we 
don't have at all "excluded" region. Moreover, it appears that it is general situation when 
we start not from GUT scale but from Planck scale (that is more "scientific"). Even if 
we work within SU(3) ®SU(2) ® U(l) MSSM and impose boundary conditions at Planck 
scale we find that in this case also the lightest sparticle is neutralino. For instance, for 
mo = 0,mi = 500Gev,A = 0,sign(fi) = - we find that in "PLANCK" ("GUT") model 
the masses are: 

m(u L ) = 1242(1118) Gev, 
m(u R ) = 1199(1077) Gev, 
m (b R ) = 1177(1072) Gev, 
m{t2) = 1120(1027) Gev, 
m {e L = 475(353) Gev, 
m(u~ L ) = 471(348) Gev, 
m(f R ) = 364(196) Gev, 
m( X2 ) = 441(401) Gev, 
m ( X4 ) = 993(890) Gev, 

m(xt) = 993(890) Gev, 
m(H°) = 1255(1102) Gev, 
m(H + ) = 1254(1101) Gev. 
We think that irrespective of validness of our concrete SU(5) model for the description 
of physics between Planck and GUT scales for small m the effects of the evolution from 
Planck to GUT scales are very important. Note also that for "Planck" model we have 



m(g) = 


= 1377(1278) Gev, 


m(d L ) 


= 1244(1120) Gev, 


m(b L ) 


= 1115(1019) Gev, 


m(tl = 


= 918(841) Gev, 


m{v L ) 


= 471(348) Gev, 


m(e R ) 


= 404(196) Gev, 


m(f L ) 


= 475(352) Gev, 


m(xi) 


= 226(205) Gev, 


m(xa) 


= 989(886) Gev, 


m(xt) 


= 441(401) Gev , 


m(h) = 


= 108(106) Gev, 


m(A) 


= 1252(1098) Gev, 



rather big splitting among the sleptons of the first two generations and the third one (in 
our program sparticle masses of the first and second generations coincide). Second lesson 
from our program is that for the case of strong violation of gaugino mass universality at 
GUT scale we can have the situation when the second neutralino and the first chargino 
being mainly SU(2) gaugino are heavier than the gluino. For instance, for tuq = 800Gev, 
mi = 200 Gev, A = 0, tan/? = 2, 5 = — 1, sign(fi) = — in "GUT" model gluino mass is 
m(g) = 599 Gev, whereas the neutralino and chargino made basically from gaugino with 
small mixing of higgsino have the following masses: 

m(xi) = 168 Gev, m(x°) = 664 Gev, m(xi) = 664 Gev. All squarks masses here are 
heavier than the gluino mass. It means that gluino can decay only to quark- ant iquark 
and LSP so we shall not have multilepton signatures and the single SUSY signature at 
LHC and TEVATRON will be multijets with missing Et- Moreover for the case of the 
model with R-parity violation of the uds type LSP will decay into two quarks that makes 
the observation of such particular scenario very difficult at LHC. 



A brief overview of the code 

A code is written on FORTRAN 77; we might not be keeping to all the standards, but 
at least, it runs OK under UNIX on CERN clusters. Double precision is used throughout 
all the calculations. Systems of algebraic and differential equations are solved by calls to 
CERNLIB |l3j subroutines DSNLEQ and DDEQMR. For neutralino matrix eigenvalues 



subroutine RSM is used (taken from |14]| ). Some checks are performed during the cal- 
culations, so for some "bad" initial values, the code skips further calculations at certain 
moment and returns with non-zero error flag (see below). Still, no checks on any physical 
meaning is made, so it is up to user to analyse output masses and error flag. Essentially 
procedure consists of solving the systems of differential equations with some of the vari- 
ables (h t ,h T ) being fixed in fact at the final scale. To find their values at the starting 
scale ( Mgut or Mpi anc k) we start with the solving the system of equations, decoupled 
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indeed for most (but not for all!) reasonable input parameters : 

\electroweak scale fixed /if, 
h T {jj) \electroweak scale fixed fl T 

in respect to x,y (values of h t ,h T at starting scale). Their values at electroweak scale 
are fixed by the values of m t ,m T ,a s . In the process of solution, we find also the scale, 
taken as \f m(ti)m(t 2 ) to change from two- to one-loop RG equations. To get the values of 
m(ti),m(t2) we solve the equation to find electroweak minimum fi. The rest of program 
flow is rather trivial arithmetics except for masses "tuning" (their values are taken not at 
common scale but each one is calculated at the scale equal to this very mass). It requires 
again the solution of equation for the value of the mass, now with the final scale taken as 
variable. 

Masses calculation is performed by the call to steering subroutine SUSYM with the calling 
sequence : 

call susym 

$(iway,AMO,AM_HALF,A_t ,tanb, AMT, csign_mu, isec,dnu, idebug, 
$amass , ierr) 

User must declare 

integer iway(3) 
real amass (42) 
character*l csign_mu 

Input parameters : 

• — iway(l) = 1 : to run from Planck to weak scale 

- iway(l) = 2 : from GUT to weak scale 

- iway(2) = 1 : normal run ( all calculations of all the masses) 

- iway(2) = 2 : Stop after h t) h b tuning. Useful in looking for possible pairs of 

(tan(P),M t ) 
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— iway(3) is not used for now 

• AMO : M at initial scale (real) 

• AM_HALF : M 1/2 at initial scale (real) 

• A_t : A t (real) 

• tanb : tan(/3) (real) 

• AMT : M t) mass of top-quark (real) 

• isec : second loop coefficient (integer); set it normally to f 

• dnu : non-universality coefficient (real) 

• idebug : produce some output if > (integer) 
Output parameters : 

• array amass : values/masses at the weak scale 

1. 9l 

2. gl 

3. gl 

4. M 1/2 (l) 

5. M 1/2 (2) 

6. M 1/2 (3) 

7. ht 

8. ht 

9. h T 
fO. A t 
11. A b 
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12. A T 

13. A 

14. A 

15. A 

16. A 

17. £>, 

18. e/ 

19. e r 

20. A 

21. T) 

22. B r 

23. T r 

24. v Tl 

25. f, 

26. f r 

27. // 

28. Ti 

29. f 2 

30. #i 

31. # 2 

32. x? 

33. X °2 

34. xl 

35. x° 



36. 


xt 


37. 


xt 


38. 


H axial 


39. 


H charged 


40. 


Hlight 


41. 


Hheavy 


42. 


Gluino 



• ierr = 

— OK 

— 1 error in input parameters 

— 10 bad fitting of h t , hb 

— 100 negative m 2 for some of sparticles 

— 1000 first approximation for fi 2 is negative 

— 2000 low and high limits for /i 2 are 0. both 

— 3000 no good solution for fi found 

— 10000 failed to find eigenvalues for neutralinos 

— 100000 failed to find eigenvalues for charginos 

To conclude, we have written "PLANCKSUSY" program for SUSY masses calcula- 
tions. The main peculiarity of our program for instance in comparison with ISASUSY 
program |T2j is that we can start from both Planck and GUT scales. We also introduced 
the non-universality parameter for gaugino masses. It appears that for very interesting 
region in the parameter space when mo is much smaller than mi if we start from the 
Planck scale we find that in this region the LSP is always neutralino unlike to the case 
when we start from GUT scale and LSP is the right-handed tau-slepton. For the case 
when we have strong deviation from gaugino mass universality it is possible to realize 
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the situation when the second neutralino and the first chargino are heavier than gluino. 
For such scenario lepton signatures for SUSY search at LHC and TEVATRON are absent 
and the single SUSY signature are multijet events with missing Et- Moreover, for the 
models with R-parity violation of uds type we shall have only jet events in final states that 
makes SUSY discovery for such scenario rather difficult. The program can be received by 
e mail:Vsevolod.Popov@cern.ch 

We thank CERN CMS Department for the hospitality during our stay at CERN where 
this paper has been finished. We are indebted to the participants of the Daniel Denegry 
CMS physics simulation seminar and espessially to Daniel for vivid discussions. We also 
indebted to I.Semenjuk for the help and consultations. 
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